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FOREWORD 


Bell Helicopter Textron Inc. (BHTI) has been conducting a study of finite element modeling of 
helicopter airframes to predict vibration. This work is being performed under U.S. Government 
Contract NAS1-17496. The contract is monitored by the NASA Langley Research Center, Structures 
Directorate. 

This report summarizes the procedure used at BHTI for predicting coupled rotor/fuselage 
vibrations with an application to the AH-1G two-bladed rotorcraft including comparisons with 
flight test vibrations. Key NASA and BHTI personnel are listed below: 


NASA Lanalev 

Panice H. Clark, Contracting Officer 

Joseph W. Owens, Contract Specialist 

John H. Cline, Technical Representative 

Raymond G. Kvaternik, Leader, Rotorcraft 

Structural Dynamics Group 


Bell Helicopter Textron Inc. 

W. Young, Manager, Research 

J. D. Cronkhite, Group Engineer, 

Structural Dynamics 

R. V. Dompka, Senior Research Engineer 
G. Sadler, Chief, Rotor Dynamics 
J. Corrigan, Rotor Dynamics Group Engineer 

J. Rogers, Rotor Dynamics Engineer 

K. S. Perry, Rotor Dynamics Engineer 
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INTRODUCTION 


The NASA Langley Research Center is sponsoring a rotorcraft structural dynamics program with the 
overall objective to establish in the United States a superior capability to utilize finite 
element analysis models for calculations to support industrial design of helicopter airframe 
structures. Viewed as a whole, the program is planned to include efforts by NASA, universities, 
and the U.S. helicopter industry. In the initial phase of the program, teams from the major U.S. 
manufacturers of helicopter airframes will apply extant finite element analysis methods to 
calculate static internal loads and vibrations of helicopter airframes of both metal and com- 
posite construction, conduct laboratory measurements of the structural behavior of these air- 
frames, and perform correlations between analysis and measurements to build up a basis upon which 
to evaluate the results of the applications. To maintain the necessary scientific observation 
and control, emphasis throughout these activities will be on advance planning, documentation of 
methods and procedures, and thorough discussion of results and experiences, all with industry- 
wide critique to allow maximum technology transfer between companies. The finite element models 
formed in this phase will then serve as the basis for the development, application, and 
evaluation of both improved modeling techniques and advanced analytical and computational 
techniques, all aimed at strengthening and enhancing the technology base which supports 
industrial design of helicopter airframe structures. Here again, procedures for mutual critique 
have been established, and these procedures call for a thorough discussion among the program 
participants of each method prior to the applications and of the results and experiences after 
the applications. The aforementioned rotorcraft structural dynamics program has been given the 
acronym DAMVIBS (Design Analysis Methods for VIBrationS). 

Under the DAMVIBS program, the four industry participants (BHTI, Boeing Helicopters, McDonnell- 
Douglas Helicopter, and Sikorsky Aircraft) are to apply existing company methods for coupled 
rotor-fuselage analysis to calculate vibrations of the AH-1G helicopter and to correlate with 
data available from an Operational Load Survey (OLS) flight test program (References 1 and 2). 
In support of this common activity, BHTI, the manufacturer of the subject aircraft, was tasked to 
prepare and provide to the other participants the data needed to independently make these 
analyses and correlations. Specifically, BHTI was tasked to: 

1. Present a detailed description of the modeling rationale and techniques used to develop 

the AH-1G NASTRAN fuselage finite element vibration model under a previous contract 

(Reference 3). A NASTRAN data deck of this model was provided to the other 
participating manufacturers. 

2. Present a detailed description of all previous correlation work used to verify the 

finite element model (two versions - stick and built-up tailboom), including the 

following: 
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a. Ground vibration tests (GVT), static deflection tests and in-flight excitation 
simulation (References 4 and 5). 

b. Application of the built-up tailboom model predictions to the static and ground 
vibration tests of Reference 4. 

c. Correlation of both models with other prior AH-1G results contained in References 
6 and 7. 

3. Describe the OLS flight test program on the AH-1G and assemble the vibration data to be 
used in the correlations. 

4. Present the AH-1G rotor system mechanical and aerodynamic coefficient data to all 
participants. 

References 1-7 were used to develop the necessary background on the FEM and flight loads data for 
the current rotor/fuselage coupling analysis task as summarized in References 8 and 9. 

This report describes work conducted by BHTI to evaluate the adequacy of current theoretical 
methods for predicting coupled rotor/fuselage vibration. The analysis methods described herein 
represent BHTI's advanced analysis rotorcraft flight simulation computer program C81 (Reference 
10) and the industry's primary structural analysis computer program NASTRAN. These analytical 
methods form the basis of an approach to helicopter dynamic analysis that has been used 
successfully at BHTI for many years. This report describes the analytical formulation of rotor 
dynamic equations, fuselage dynamic equations, coupling between the rotor and fuselage, and 
solutions to the total system. The coupled analysis is applied to an AH-1G two-bladed rotor 

system and results compared with measured OLS flight test vibrations. 
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2. DESCRIPTION OF ROTOR/FUSELAGE 

COUPLING METHOD 



GENERAL PROCEDURE 


The general procedure used in this study for calculation of AH-1G vibration characteristics is 
depicted on the following page. A Myklestad-type of analysis is used to calculate rotating 
elastic blade modal properties. A NASTRAN finite element method is used to form isolated 
fuselage modal characteristics. C81 is then used to couple the rotor and airframe math models, 
and then calculate rotor hub loads. Finally, the calculated hub loads are used as a forcing 
function on the NASTRAN finite element model to calculate the vibration responses. 

The Myklestad family of programs has been used at BHTI to calculate helicopter rotor blade 
natural frequencies and mode shapes for many years. The capabilities of the program have been 
modified to include a complete elastic and inertial representation of the blade, pitch control 
systems, and pylon impedances. Recent modifications have been made to enhance user convenience, 
to give more accurate results, and to model more advanced rotor configurations. 

C81 is a comprehensive rotorcraft flight simulation program used to calculate the aeroelastic 
response of the coupled rotor/fuselage system. The structural analysis is based on a modal 
technique while rotor aerodynamics are modeled using strip theory and bivariant tables to 
represent stall and compressibility effects over the rotor disk. 
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GENERAL PROCEDURE 












C81 HISTORICAL DEVELOPMENT 


Initially, C81 rotor performance and handling qualities were based on an actuator disk 
representation with a simple 6 degree-of-freedom rigid fuselage. In 1968, rotor aerodynamics 
graduated to strip theory representations with analytic functions dependent on blade radius, 
location, and time representing air loads. In addition, the fuselage and rotor were dynamically 
coupled. In 1971, Cl and Cg coefficients were added to the strip theory as well as an 
aeroelastic rotor model and rigid pylon representation. 1973 saw the first acceptance of a C81 
deck in an analytic competition for rotorcraft design. In 1977, Floquet theory was added to 
address the stability of systems whose motion was governed by differential equations with 
periodic coefficients (e.g. ground resonance effects). Optimization capabilities were added in 
1983 to improve the usefulness of C81 in design and in 1985 the capability to handle 10 elastic 
fuselage modes (in addition to six rigid modes) was included in the coupled rotor-fuselage 
program. 
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C81 HISTORICAL DEVELOPMENT 


1965 ROTOR AERODYNAMIC MODEL 

1968 COUPLED ROTOR /FUSELAGE (RIGID) 

1971 AEROELASTIC ROTOR 

1973 REQUIRED IN AAH/UTTAS COMPETITION 

1976 DISSIMILAR BLADES 

1977 FLOQUET THEORY ANALYSIS 

1983 OPTIMIZATION 

1985 ELASTIC FUSELAGE EFFECTS 
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C81 ANALYSIS CAPABILITY 


C81 has the capability to model a wide variety of aircraft and hub configurations as shown on the 
following page. BHTI's comprehensive rotorcraft dynamics analysis code C81 is capable of 
modeling the following components of a rotorcraft: A fuselage; two rotors, each with a modal 
pylon, aeroelastic blades, and a nacelle; a wing; four stabilizing surfaces, none of which must 
be purely vertical or horizontal; four external stores or aerodynamic brakes; a nonlinear, 
coupled control system including a collective bobweight, stability and control augmentation 
system, and maneuver autopilot simulator; two jets; and a weapon. 
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C81 ANALYSIS CAPABILITY 
AIRCRAFT CONFIGURATIONS 


• CONVENTIONAL MR/TR • TILT ROTOR 

• TANDEM • COAXIAL 

• SIDE BY SIDE • COMPOUND 

• AIRPLANE • DIRIGIBLE 

• WINDMILL • WIND TUNNEL 


HUB CONFIGURATIONS 

• GIMBALLED 

• TEETERING 

• ARTICULATED 

• HINGELESS 
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C81 COORDINATE SYSTEMS 


A C81 analysis employs five coordinate systems that describe the behavior of various rotorcraft 
components. Each coordinate system is denoted by the following subscript notation. 

Subscript 

Hone 1. Ground Reference - A non-rotating coordinate system taken to be the inertial 

reference system. 

f 2. Fuselage Reference - A non-rotating coordinate system centered at the fuselage 

center of gravity. 

m 3 . Mast Reference - A non-rotating coordinate system centered at the top of the 

mast. 

h 4 . Hub Reference - A rotating coordinate system that shares the same origin as the 

mast reference system. 

b 5. Blade Reference - The origin of the blade rotating coordinate system is at the 

inboard end of the feathering bearings. 
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BHTI ELASTIC ROTOR BLADE DYNAMIC ANALYSIS 


Elastic rotor blade dynamic characteristics are calculated for use in C81 by a Myklestad - based 
analysis implemented in a computer code called DNAM06. Blade modal natural frequencies, 
generalized inertias, and mode shapes are calculated by DNAM06 in rotating coordinates. All 
linear mass and spring terms are included in the DNAM06 analysis. Terms that were deleted from 
the DNAM06 analysis in order to simplify the Myklestad analysis, such as the Coriolis 
acceleration terms, are included in the forcing function of the C81 rotor analysis. Nonlinear 
terms such as flapping springs and flap-lag-torsional moments are also handled by the C81 
analysis, as are the coupling effects of the rotor hub's motion. 

The following figure lists the assumptions which underlie the linear blade equations of motion in 
the Myklestad analysis. 
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BHTI ELASTIC ROTOR BLADE DYNAMIC ANALYSIS 


ASSUMPTIONS: 

• BLADE CROSS SECTIONS ARE NOT SYMMETRIC. 

• ELASTIC AXIS IS A PIECEWISE PARALLEL STRAIGHT LINE. 

• STRUCTURAL PRINCIPAL AXES AND MASS PRINCIPAL AXES ARE PARALLEL, BUT 
NOT COINCIDENT. 

• EFFECTS OF LINEAR AND ANGULAR DISPLACEMENTS AND ACCELERATIONS OF 
THE ROTOR HUB ARE INCLUDED. 

• EFFECTS OF STEADY BLADE FEATHERING MOTION ARE INCLUDED. 

• SHEAR DEFORMATION IS NEGLECTED. 

• PRECONE AND UNDERSLING ARE NOT INCLUDED. 
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BHTI ELASTIC ROTOR BLADE REPRESENTATION 


The DNAMQ6 analysis uses a finite element transfer matrix approach to represent blade properties. 
The rotor blade is represented by a series of lumped, rigid 3-D dumbell inertias connected by 
untwisted, massless, elastic beams. Built-in twist of the blade is introduced incrementally at 
the blade stations where the inertias are located. Rotating fully-coupled inplane, out-of-plane, 
and torsional deflections of the blade are also considered. Radial extension of the blade is 
neglected, however. 

A representative blade element is shown in the following figure. 
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BHTI ELASTIC ROTOR BLADE REPRESENTATION 


3-D INERTIA DUMBELL 
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MYKLESTAD STATE VECTOR SOLUTION 


DNAM06 uses the state vector and transfer matrix type of formulation developed by Myklestad. The 
state vectors consist of two linear displacements, three angular rotations, and the shears and 
moments corresponding to these displacements and rotations. 

The state vector contains the following quantities: 


u 

i 

i 

w 

i 


- inplane displacement 

- out-of-plane displacement 

- in-plane slope 

- out-of-plane slope 


V x -j - in-plane shear 
V 2 i - out-of-plane shear 
M z j - out-of-plane moment 
M x i - in-plane moment 

- torsional deflection about y axis 
Myi - torsional moment about y axis 

The general form of the transfer matrix equation is given on the following page. 
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MYKLESTAD STATE VECTOR SOLUTION 


{SJ = 


WHERE: 

{Sj} - STATE VECTOR AT STATION I 

{Mj} - TRANSFER MATRIX ACROSS INERTIA DUMBELL 

[fj] - TRANSFER MATRIX ACROSS ELASTIC ELEMENT OF LENGTH L f 
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ROTOR DYNAMIC ANALYSIS BY THE MODAL TECHNIQUE 


The time variant aeroelastic rotor representation in C81 is based on a modal analysis approach. 
Some of the assumptions contained in this analysis are presented on the following page. 
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ROTOR DYNAMIC ANALYSIS BY THE MODAL TECHNIQUE 


ASSUMPTIONS: 

• THE ROTOR BLADE IS DIVIDED INTO THE SAME RADIAL SEGMENTS FOR BOTH 
AERODYNAMIC AND DYNAMIC CALCULATIONS. 

• EACH SEGMENT FACE HAS THREE DEGREES OF FREEDOM THAT ARE USED IN 
THE GENERALIZED FORCE CALCULATION, (w - OUT OF PLANE, u - IN PLANE, <t> - 
TWIST ABOUT Y AXIS). 

• DNAM06 (ELASTIC BLADE ANALYSIS) SUPPLIES THE NORMAL MODES THAT 
DESCRIBE u, w, AND 4> DISPLACEMENTS FOR EACH SEGMENT FACE OF EVERY 
BLADE MODE. 

• LINEAR INTERPOLATION IS USED TO DEFINE MODE SHAPES BETWEEN TWO 
ADJACENT FACES. 

• EFFECTS OF STEADY BLADE FEATHERING MOTION ARE INCLUDED. 

• 20 SEGMENTS PER BLADE, MAXIMUM. 

• 11 INPUT BLADE MODES (FROM DNAM06) PER ROTOR, MAXIMUM. 

• 2 ROTORS, MAXIMUM. 

• 7 BLADES PER ROTOR, MAXIMUM. 
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ROTOR DYNAMIC ANALYSIS BY THE MODAL TECHNIQUE (Concluded) 


The C81 rotor analysis is designed to handle the fully-coupled blade mode shapes calculated by 
DNAM06. These mode shapes have several attributes that are important in the modal technique. 
Each mode shape is a solution to the coupled differential equations for the free vibration of the 
total rotor system. The solution, or mode shape, is obtained by deletion of all velocity 
dependent and nonlinear terms. The collection of mode shapes forms an orthogonal set. Each mode 
shape has a natural frequency, and it satisfies the appropriate set of boundary conditions 
imposed on the blade. 
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ROTOR DYNAMIC ANALYSIS BY THE MODAL 
TECHNIQUE (CONCLUDED) 


ATTRIBUTES OF BLADE MODE SHAPES: 

• EACH IS A SOLUTION TO THE LINEAR PORTION OF THE 
COUPLED DIFFERENTIAL EQUATIONS OF FREE VIBRATION 

• EACH HAS A NATURAL FREQUENCY 

• THE COLLECTION FORMS AN ORTHOGONAL SET 

• EACH SATISFIES IMPOSED BOUNDARY CONDITIONS 
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ROTOR BLADE ROOT BOUNDARY CONDITIONS 


The behavior of the rotor hub can be related to the blade mode shapes that are used to describe 
the entire rotor system by the boundary conditions presented on the next page. A hub impedance 
model is used to include the effects of an isotropic support system in the blade modes. For 
collective modes, out-of-plane motion is restrained by one input spring constant and one input 
mass. In-plane slope changes are opposed by a torsional spring. For cyclic modes, the in-plane 
motion is restrained by one input spring and one input mass. 
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ROTOR BLADE ROOT BOUNDARY CONDITIONS 


MOTION TWO BLADED 

ROTOR 


MODE TYPE 

OUT-OF-PLANE 

IN PLANE 

TORSION 

RESPONSE 

COLLECTIVE 

CANTILEVER 

PINNED 

CANTILEVER 

0, 2, 4, 6p 

CYCLIC 

PINNED 

CANTILEVER 

CANTILEVER 

1,3, 5p 
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SOLUTION FORM OF ROTOR DYNAMIC SYSTEM 


A separation of variables approach is employed in the solution of the coupled equations as 
depicted in the figure. The independent variables are blade radial location y, and time t. The 
deflections u n , w n , and 4> n of the n th elastic blade mode shape are only a function of blade 
radial position y. On the left hand side of the equation, u, w, and <j> are the total elastic 
deformation of the blade which are dependent upon both blade radial position and time. The modal 
participation factor, 8 n , for the n^* 1 mode is only a function of time. 


26 



SOLUTION FORM OF ROTOR DYNAMIC SYSTEM 


SEPARATION OF VARIABLES 


U (y, t) | 

NM 

u „ ( - v) 

w(y,t) 

> — V 

j W n ( >' ) 

<t> (y, D 

n - 1 

4> n ( y> 


5 (t) 

n 


u w <j> - COMPONENTS OF NTH BLADE MODE SHAPE. 

n* n Mi 

«„ - PARTICIPATION FACTOR FOR NTH BLADE MODE. 

TOTAL BLADE ELASTIC DEFORMATION. 
nm- NUMBER OF BLADE MODES. 


27 



ROTOR BLADE AERODYNAMIC REPRESENTATION 


A brief summary of C81 rotor aerodynamic capabilities is listed below: 

1. Maximum of 20 radial aerodynamic segments. 

2. The airfoil sectional C-j, C d , and C,,, are tabular functions of Mach number and angle of 
attack. 

3. Momentum theory induced velocity is calculated based either on equations internal to the 
program or user input tables. 

4. Unsteady aerodynamic options include Theodorsen and Carta theories. 

5. The effects of blade elasticity are included. 

The rotor blade aerodynamic reference system is shown below. 
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MODAL FORM OF ROTOR AEROELASTIC EQUATIONS 


Application of the separation of variables technique to the solution of the rotor blade equations 
of motion results in the modal blade equations presented in the figure. Here, F u is the applied 
inplane force, F^ is the applied out-of-plane force, M<t, is the applied twisting moment and I n is 
the generalized inertia of the n* h blade mode. F n is the forcing function due to the aerodynamic 
forces and inertial forces, including those forces which were deleted from the Myklestad 
analysis. It should be noted that the left hand side of the rotor equations is uncoupled due to 
the orthogonality of the mode shapes obtained from DNAM06. 
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C81 ROTOR AERODYNAMIC MODULE 


The forcing function terms A,,, A w , 
motion are calculated in the l 81 ro 


on the right 


rotor aerodynamic module. 


side of the rotor blade equations 
Representative inputs and outputs 


this module are shown in the diagram below. 
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C81 ROTOR AERODYNAMIC MODULE 
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GENERAL FORM OF BLADE EQUATIONS 


The set of previously presented modal equations is sufficient to describe the time variant aero- 
elastic response of the rotor blade because of two reasons: 

1. The linear effects of blade mass, blade elasticity, and blade geometry are inherently 
included in the calculated blade mode shapes and natural frequencies. 

2. Aerodynamic and aeroelastic effects are included in the forcing terms on the right hand side 

of the blade modal equation. Dynamic terms which were deleted from the blade natural fre- 
quency and mode shape calculation are also included on the right hand side of this equation. 

These dynamic terms arise from angular velocities and accelerations . of the rotating 

coordinate system, from linear rotor hub accelerations, and from inertial terms such as 
Coriolis accelerations. 

The following figure presents the general form of the blade equation of motion along with an 
indication of the analysis in which each of the terms is calculated. 
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GENERAL FORM OF BLADE EQUATIONS 



CALCULATED IN MYKLESTAD CALCULATED IN C81 

BLADE MODE ANALYSIS 
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FUSELAGE DYNAMIC ANALYSIS 


The aeroelastic behavior of a rotor is highly dependent upon the dynamic characteristics of the 
fuselage to which it is attached. The fuselage (or nonrotating components) is represented in C81 
by the previously discussed modal technique. This modal form of the fuselage equations, was 
selected in order to provide general motions at the rotor hub with a minimum number of additional 
equations (thus reducing computational requirements). This technique allows the use of fuselage 
modes based on extremely complex and detailed fuselage/support system models which can be 
obtained from structural analysis finite element codes such as NASTRAN. Each fuselage mode shape 
can have three linear displacements and three angular rotations at the rotor hub as part of the 
modal information. The modal information also contains fuselage natural frequencies, generalized 
inertias, and damping where applicable. The C81 analysis is currently designed to handle up to 
ten elastic fuselage modal equations. For convenience, C81 has an option to input fuselage modes 
that were calculated with or without full rotor mass (included as a point mass), for use in 
coupled rotor/fuselage dynamic analysis. A simplifying assumption used by the C81 fuselage 
dynamic analysis is that the forcing function comes from the rotor, and that there are no other 
aerodynamic or inertial loads applied to the fuselage model. 
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FUSELAGE DYNAMIC ANALYSIS 


• NO OSCILLATORY AERODYNAMIC FORCES ARE APPLIED TO AERODYNAMIC 
SURFACES 

• HUB DEGREES OF FREEDOM AND HUB LOADS (FROM THE ROTOR) ARE 
USED IN THE GENERALIZED FORCE CALCULATIONS. 

(THREE LINEAR DISPLACEMENTS: X j# Y jf Z ] AND THREE 

ANGULAR ROTATIONS: e xj , e yjf e zj OF THE HUB) 

• NASTRAN (FINITE ELEMENT ELASTIC FUSELAGE ANALYSIS) SUPPLIES THE 
NORMAL MODE DATA 

• 10 INPUT ELASTIC FUSELAGE MODES (FROM NASTRAN), MAXIMUM 





FUSELAGE MODAL EQUATIONS OF MOTION 
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FUSELAGE MODAL 


p + 2 £ . o 
j J 




FUSELAGE MODAL FORCING FUNCTION 


The modal forcing function F; on the right hand side of the fuselage modal equation of motion is 
given below. The quantities J V mx , V m w, V m2 , M mx , M m y, and M mz are the shear and moment components 
at the top of the rotor mast. The^e shears and moments are calculated by a modal displacement 
technique, which combines the hub shear and moment coefficients obtained from the Myklestad 
analysis and the solution of the rotor's modal equations of motion. Fmx. F m y. and F mz are force 
components due to any translation of the rotor mass not included in either tne rotor or fuselage 
analysis. Since there is no radial degree of freedom in the rotor dynamics model, Fmx and F m y 
also include the corrections needed to account for the inertial forces associated with radial 


foreshortening. 
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FUSELAGE MODAL FORCING FUNCTION 


F = X V + Y • V + Z • V 

1 j ^mj mx 1 mj w my ‘“mj mz 

+ ®mxj ^mx + 0 myj ^my + ®mzj ^mz 
^mj ^mx * ^mj ^my * ^mj ^mz 


OR, 


= F 0 j ^mj ^mx + ^mj ^my ^mj ^mz 
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TRANSFORMATION OF FUSELAGE COORDINATES 


The fuselage mode shape is defined in terms of the following displacements at the top of the 
mast expressed with respect to the fuselage reference system: 

xj - in Xf direction 
yj - in yf direction 
Zj - in Zf direction 
e X j - about Xf axis 
0 y j - about yf axis 
0 2 j - about Zf axis 

In order to perform the coupled rotor/fuselage analysis, the fuselage mode shapes must be 
transformed into the mast reference system. The transformation equation is presented below, 
where the subscript "m" refers to the mast reference system and no subscript refers to the 
fuselage reference system. This transformation references the fuselage equations to the same 
coordinate system as the rotor equations. 
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TRANSFORMATION OF FUSELAGE COORDINATES 


(X ro j» ^mj* ^mj) = ^j' ^j' 

( e xmj' e ymj» e zmj) = D"m/fl ^ e xj' G yj' 

WHERE: [T m/f ] = FUSELAGE-TO-MAST COORDINATE TRANSFORMATION MATRIX. 


43 



ROTOR HUB ACCELERATION COMPONENTS 


The total linear acceleration at the rotor hub is thejsum of fuselage elastic contributions and 
fuselage rigid body contributions. The terms x* m , y m , and z m are the rigid body fuselage 
acceleration components written with respect to the mast reference axis system. The terms axm« 

.2 -ft . 

aym»azm represent the hub accelerations which come from the term w ___ on pg 35. NP is the 
number of fuselage elastic modes used in the analysis. 
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ROTOR HUB ACCELERATION COMPONENTS 


(a 


xrrv ^ 



NP .. 

^ P j (*mj» Yrrij' ^mj ) + (^nv Ym'^m) 

] = 1 


• NP = NUMBER OF FUSELAGE ELASTIC MODES USED IN ANALYSIS 

• MAST REFERENCED (NON -ROTATING) COORDINATES 

• X m# y m# Z m - RIGID FUSELAGE ACCELERATION COMPONENTS 
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IN-PLANE INERTIAL FORCES DUE TO ROTOR HUB MOTION 


The inplane rotor hub accelerations, a xm and ay m . produce in-plane forces, Fj and Fr perpendicular 
and parallel to the blade, respectively. The inertial loading is separated in this fashion because 
for an elastic rotor the Fj force is absorbed into the rotor modal equations as part of the forcing 
function for the rotor inplane degree of freedom, while the Fr term has to be handled separately 
since there is no radial degree of freedom in the blade modes. 
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IN-PLANE INERTIAL FORCES DUE TO ROTOR HUB MOTION 


| 



i 



INERTIAL HUB FORCES DUE TO ROTOR HUB MOTION 


After integrating over each blade and summing the contribution of all blades, the inplane rotor 
hub forces F ym and F^ are given by the following equations. For an elastic rotor analysis the 
Fj terms in cne equation must be omitted because their effect will be already accounted for in 


the modal analysis. 
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INERTIAL HUB FORCES DUE TO ROTOR HUB MOTION 


F = 

xm 


NB 

E 

i = 1 


/ (-dFjSIN'F. -dF R COSV.jdy 
o 


F 

yin 


NB 

S 

i = 1 


/ ^-dFjCOSV. +dF R SINM».^dy 
o 


Fi omitted for elastic rotor analysis 
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IN-PLANE INERTIA FORCES FOR AN ELASTIC ROTOR 


For an elastic rotor, the inplane force perpendicular to the rotor blade, Fj t is already included 
in the forcing function that is applied to the elastic rotor blade modal equations, and therefore 
the hub loads due to these terms are computed by summation of the load over the number of blade 
modes. The inplane force in the radial direction, Fr, has to be integrated separately and 

included as a hub force in the fuselage equations. Fr has to be handled n^vklStad 

stated earlier, there is no radial degree of freedom in the elastic blade model in Myklestad. 

The hub shears due to Fr are presented on the following page. 

It should be noticed that the out-of-plane inertial force included arising from 

included. For an elastic rotor, a distributed inertia force is included as part of _the roto 

forcing function, and the vertical shear that is applied at the rotor hub is computed from the 

modal response of the rotor. 

Also, the inertia forces on the rotor caused by angular velocities and angular accelerations of 
the fuselage are also included in the rotor forcing function. These effects are transmitted to 
the fuselage at the rotor hub by means of the net beamwise bending moments computed from the 

modal response of the rotor. 
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WHERE: 


IN-PLANE INERTIA FORCES FOR AN ELASTIC ROTOR 


NB NP 

P = -M b 2 E p(x COSV-y •COSH'SIN'P) 

xm a *J m) » 1 1 

i = i j = i 
NB 

- M b e <x m COS 2 V. - y m COS •P.SIN V.) 
i = 1 

NB NP 

F = M « 2 S p.(x . COS V. SINV. — y .SIN 2 V.) 

ym B J n»J i i mj i 

»=1 j=l 

NB 

+ m b 2 (x m COS V. SIN 4». - y m SIN 2 V.) 
i= 1 

M b = MASS OF ONE BLADE. 

NB = NUMBER OF BLADES 
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MODAL FORM OF FUSELAGE EQUATIONS OF MOTION 


The fuselage modal equations of motion, which include rotor mass, are given on the following page 
for a two-bladed elastic rotor. The terms with a single underline are needed to account for the 
rotor mass that was included in the NASTRAN analysis. The double underlined terms are the result 
of the rotor's radial inertial terms which are not handled by the rotor modal analysis since 
there is no radial degree of freedom. The terms on the right hand side of the equation are due 
to the rigid body motions (%, ym* zm) of the airframe. 
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MODAL FORM OF FUSELAGE EQUATIONS OF MOTION 

• TWO BLADED ROTOR 

• ROTOR MASS INCLUDED 


2M, 


NP 


’VST 


s Pi 
J i = 1 


x .* . (COS 2 ¥ - 1 ) + y . y . (SIN 2 V - 1 ) 

mi mj nil * mj _ 


_Z z - COS ¥ SIN 'P (y x +x y .) 

mi tnj ,m tn J mi m J 


+ 2 S“ J p y 


F . 

. ~ 2 _ °J 

+ p — 1,1 


2M„ 


J GI GI 
j J 


i x . (COS* V - 1) + y m y (SIN 2 W - l) 

m mj w mj 


-i z -COS'I'SINH'(yx + x y ) 

m mj m tnj m mj 


WHERE: Mb = MASS OF ONE BLADE 
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SOLUTION SCHEME FOR DYNAMICALLY COUPLED EQUATIONS OF MOTION 


The following sequence of calculations is used to solve the coupled rotor/fuselage dynamic sys- 
tem. This sequence was developed to provide a consistent set of equations without the need to 
solve the full set of equations simultaneously. 

1. Find aerodynamic loading for rotors, empennage, and fuselage. These forces depend only upon 
displacements and velocities. 

2. Compute that portion of the rotor modal forcing functions not dependent on fuselage accel- 
erations (such terms as nonlinear flapping springs and dampers and unsteady aerodynamic 
effects) . 

3. Solve for rigid body fuselage accelerations due to rotor excitation. 

4. Solve for the accelerations of the fuselage generalized coordinates due to rotor excitation. 

5. Add the inertia loads caused by fuselage motion to the rotor modal forcing function. 

6. Solve for rotor accelerations. 

Hammings Predictor-Corrector Method of numerical integration is used to integrate the equations 
of motion. 
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STOP 


SOLUTION SCHEME FOR DYNAMICALLY COUPLED 

EQUATIONS OF MOTION 
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9S 




APPLICATION TO THE AH-1G OPERATIONAL LOAD 
SURVEY (OLS) FLIGHT TEST PROGRAM 
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Summary of Required C81 Inputs 


The inputs to C81 that were necessary to calculate AH-1G OLS hub loads are presented in the figure on the 
opposite page. These inputs consist of rotor blade modal data, fuselage modal data, and basic C81 input 
data. The rotor blade modal data is calculated by the Myklestad computer program which used OLS inputs 
from Reference 2 to describe the physical rotor blade properties. Blade modal data calculated by 
Myklestad and used as inputs to C81 consists of natural frequencies, generalized inertias, and mode 
shapes. Isolated fuselage modal data is calculated by NASTRAN (based on OLS data) and these data (in the 
form of natural frequencies, generalized inertias, and mode shapes) are used in C81 to describe fuselage 
dynamic behavior. All other basic C81 inputs were taken from Reference 8 and these inputs, combined with 
the rotor and fuselage modal input data, comprise the total AH-1G OLS C81 input deck. Each of these three 
sets of data are discussed on the following pages. 
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SUMMARY OF REQUIRED C81 INPUTS 
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Rotor Blade C81 Input Data 


The rotor blade modal data required by C81 is calculated by the BHTI Myklestad (DNAM06) computer program. 
OLS data is used as input to this program to describe the physical rotor blade properties. This program 
then calculates the blade natural frequencies, generalized inertias, and mode shapes for input to C81. 
For the present analysis, nine modes were chosen to represent the flexible blade characteristics. Five of 
these modes are of the cyclic variety which have pinned out of plane, cantilever inplane, and cantilever 
torsional blade boundary conditions. The other four modes are collective in nature and have cantilevered 
out-of-plane, pinned inplane, and cantilevered torsional blade boundary conditions. 
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ROTOR BLADE C81 INPUT DATA 


CALCULATED BY MYKLESTAD (DNAM06) 
MODE 

1ST OUT-OF-PLANE BENDING 
1ST IN-PLANE BENDING 
1ST TORSION 

2ND OUT-OF-PLANE BENDING 
3RD OUT-OF-PLANE BENDING 

1ST OUT-OF-PLANE BENDING 
1ST TORSION 

2ND OUT-OF-PLANE BENDING 
3RD OUT-OF-PLANE BENDING 

PLUS RESPECTIVE GENERALIZED 


FREQUENCY (p) TYPE 


1.0000 CYCLIC 

1.3024 CYCLIC 

2.3781 CYCLIC 

2.7489 CYCLIC 

4.5301 CYCLIC 


COLLECTIVE 

COLLECTIVE 

COLLECTIVE 

COLLECTIVE 


1.0423 
2.3356 
2.9016 
4.7384 

INERTIAS AND MODE SHAPE DATA 



Fuselage C81 Input Data 


NASTRAN was used to calculate the fuselage modal data for use in C81. The full rotor weight, which 
includes the rotor, hub, and R-MUX (OLS rotating mu1tip 1ex T ^ nstru ™®^Vd° NASTRAN then expressed 

^ formations i sten t^w Ulf cT i npSf r^frl^^ 

transformed ,nto the Ml coordinate' systeS. C81 is capable of handling ten elastic 

thin in addition to six rlqid body inodes of the fuselage calculated by NASTRAN, ten elastic mooes were 

used* to represent the fuselage in the C81 analysis. These modes _ are th JnTtte 
that are most important to a rotor dynamics analysis, namely the pylon pitch and roll modes, 

fuselage first bending and torsion modes. 
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FUSELAGE C81 INPUT DATA 


• CALCULATED BY NASTRAN 


MODE 

FREOUENCY (HZ) 

PYLON PITCH 

2.99 

PYLON ROLL 

3.86 

1ST FUSELAGE LATERAL BENDING 

7.12 

1ST FUSELAGE VERTICAL BENDING 

7.96 

1ST FUSELAGE TORSION 

16.03 

2ND FUSELAGE VERTICAL BENDING 

17.22 

2ND FUSELAGE LATERAL BENDING 

17.77 

FUSELAGE ROLL/ ENGINE LATERAL 

19.26 

MAIN ROTOR MAST LATERAL BENDING 

25.59 

MAIN ROTOR MAST FORE -AFT BENDING 

27.10 


• PLUS RESPECTIVE GENERALIZED INERTIAS AND MODE SHAPES 
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AH-1G OLS C81 TRIM OPTIONS 


Two helicopter C81 trim options were employed in this study. The first technique is called "trim to 
cyclic - fall through trim". This is basically a "wind tunnel" type of trim and therefore rigid body 
effects are not included. These rigid body dynamics are important, however, for accurate loads 

calculations. To account for rigid body effects, the helicopter is allowed to fall through trim for a 
length of time corresponding to ten rotor revolutions. Ten rotor revolutions were chosen so that no 
significant changes between measured and calculated flight conditions would develop. The key point here 
is that measured OLS control position data is used as input. The second trim procedure is called "full 
aircraft trim" and is representative of the type of analysis that is used during design phases of 
helicopter development. Here, only measured OLS flight conditions are used as inputs to C81 and C81 
calculates all control positions required to trim the helicopter. The C81 input deck was then submitted 
for execution and calculation of rotor hub loads. 
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AH-1G OLS C81 TRIM OPTIONS 



CALCULATED 
ROTOR BLADE LOADS 
ROTOR HUB LOADS 




Basic C81 Input Data 




which are found in Reference 8. 
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BASIC C81 INPUT DATA 


• DATA OBTAINED FROM REFERENCE 2 (OLS DATA) THAT IS USED FOR 
CYCLIC -FALL THROUGH TRIM" OPTION 

1) ROTOR 

2) FUSELAGE 

3) ROTOR AERODYNAMICS 

4) FLIGHT CONDITIONS 

• ADDITIONAL DATA OBTAINED FROM REFERENCE 8 THAT IS USED FOR 
AIRCRAFT TRIM" OPTION 

5) TAIL ROTOR 

6) FUSELAGE AERODYNAMICS 

7) WINGS 

8) STABILIZING SURFACES 

9) CONTROLS 


TRIM TO 


"FULL 
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AH-1G OLS TRIM CORRELATION 


Comparisons between the measured blade feathering and flapping angles and the C81 calculated feathering 
and flapping angles are presented in the following two figures, while a comparison of tjje ^asured 
fuselage pitch and the C81 computed pitch angles is given in the third figure. The i test data are 
represented by the open symbols while the C81 results are represented by the solid and dashed ‘15! 
solid line is the result of the "trim to cyclic-fall through trim" option. The dashed line is the result 
of using the "full aircraft trim" option. 
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Harmonic Analysis of Calculated Time History Data 


C81 computes all loads in time history format and thus requires a harmonic analysis to separate out the 
magnitudes of the various harmonic loads found in the system. This analysis, is "self-contained in C81 
and represented by the equation on the next page. C81 determines the coefficients ao. ak» bk where k in 
this case is 2 , 4, or 6. The term a 0 represents the steady load component and the ak and bk are the 
oscillatory cosine and sine components of the time history data. Measured OLS data are harmonically 
analyzed in a similar manner. 
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HARMONIC ANALYSIS OF CALCULATED 
TIME HISTORY DATA 


• CALCULATE 2,4,6 PER REV HUB LOADS FROM CALCULATED TIME HISTORY DATA. 


N 

fit.) = a o + y 
k = 1 


a k cos (2 II k w t) 4- fc> k sin (2 II k co t) 


• COMPARE COMPONENTS TO MEASURED OLS TEST DATA 
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AH- 16 OLS Rotor Blade Loads Correlation 


Comparisons between measured OLS test data and C81 calculated rotor blade loads are presented in the 
figures on the next three pages. Some measure of the capability of C81 to predict blade loads (and thus 
hub loads) can be gleaned from this comparison. It is these blade loads that sum together to form the 
overall hub loads that will be used in the fuselage vibration prediction. Results are presented for three 
airspeeds: 67, 114, and 142 KTAS. This covers the low, medium, and high speed flight regimes of the AH- 
1G OLS data. The rotor blade beam bending moments, chord bending moments, and torsional moments are 
presented as a function of blade radial station at each airspeed considered. The test data are 
represented by the solid symbols and the “trim to cyclic-fall through trim" and "full aircraft trim" 
analyses are represented by the solid and dashed lines, respectively. 
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AH-1G OLS ROTOR BLADE LOADS CORRELATION (CONT'D) 


A test data 

TRIM TO CYCLIC-FALL THROUGH TRIM 

FULL AIRCRAFT TRIM 



RADIAL POSITION (%) 
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C81 AH-1G OLS Hub Shear Predictions 
(Trim to Cyclic-Fall through Trim) 


The table shows the hub shears that were predicted by C81 using the "trim to cyclic-fall through trim" 
procedure outlined earlier. Here, the 2 , 4, and 6-per-rev sine and cosine components are tabulated in the 
x, y, and z directions for each of the six airspeeds considered. These shears are expressed with respect 
to the C81 mast coordinate system and thus a transformation must be performed before they can be used as a 
forcing function on a NASTRAN model. 
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C81 AH-1G OLS HUB SHEAR PREDICTIONS 
(TRIM TO CYCLIC - FALL THROUGH TRIM) 


Airspeed and 
Direction 

2p (ib) 

SINE COSINE 

4p <lb) 

SINE COSINE 

6 P (Ib) 

SINE COSINE 

67 kn 

x shear 

-76.309 

-905.036 

-124.780 

382.435 

56.418 

-75.547 

y shear 

862.029 

-210.886 

-348.996 

-125.501 

67.481 

49.302 

z shear 

429.443 

-344.259 

150.856 

-211.473 

91.159 

-31.733 

85 kn 

x shear 

-295.477 

-971.699 

-23.636 

484.145 

48.578 

-149.256 

y shear 

892.615 

-311.216 

-400.849 

5.118 

124.443 

39.693 

z shear 

648.470 

55.914 

24.282 

012.215 

36.150 

5.981 

101 kn 

x shear 

-251.855 

-1117.387 

-108.144 

552.523 

89.987 

-181.971 

y shear 

986.355 

-276.693 

-447.157 

-37.646 

150.366 

73.669 

z shear 

734.793 

471.520 

-104.211 

58.411 

6.130 

39.200 

114 kn 

x shear 

-284.161 

-1317.079 

-120.129 

658.103 

113.405 

-203.682 

y shear 

1112.491 

-266.761 

-552.942 

-57.474 

167.856 

89.103 

z shear 

797.733 

764.088 

-140.148 

9.268 

-4.255 

66.954 

128 kn 

x shear 

-111.140 

-1492.782 

-233.927 

726.946 

162.504 

-199.739 

y shear 

1196.045 

-102.522 

-628.049 

-185.606 

165.572 

130.811 

z shear 

827.044 

1075.261 

-174.716 

-132.391 

-10.406 

16.668 

142 kn 

x shear 

-1298.816 

-1383.017 

419.215 

821.964 

5.539 

-242.734 

y shear 

1065.050 

-845.804 

-706.013 

352.071 

209.111 

-19.726 

z shear 

764.981 

1487.320 

97.327 

-181.859 

-24.724 

-19.889 
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C81 AH-1G OLS Hub Shear Predictions 
(Full Aircraft Trim) 


The table shows the hub shears 
outlined earlier. Here, the 2, 4 
z directions for each of the six 
mast coordinate system and thus 
function on a NASTRAN model. 


that were predicted by C81 using the "full aircraft trim" procedure 
, and 6 per rev sine and cosine components are tabulated in the x, y, and 
airspeeds considered. These shears are expressed with respect to the C81 
a transformation must be performed before they can be used as a forcing 
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C81 AH-1G OLS HUB SHEAR PREDICTIONS 
(FULL AIRCRAFT TRIM) 


Airspeed and 
Direction 

2p (lb) 

SINE COSINE 

4p (lb) 

SINE COSINE 

6p (lb) 

SINE COSINE 

67 kn 
x shear 

-449.645 

-576.629 

172.743 

194.822 

-46.094 

-54.653 

y shear 

733.817 

-360.235 

-157.309 

128.493 

39.581 

-31.745 

z shear 

587.896 

-395.343 

259.044 

-218.556 

108.024 

59.815 

85 kn 

x shear 

-467.407 

-708.048 

124.455 

283.068 

17.687 

-100.732 

y shear 

799.933 

-350.201 

-235.761 

154.176 

74.187 

6.468 

z shear 

607.275 

196.889 

6.153 

-47.855 

14.830 

55.332 

101 kn 

x shear 

-499.956 

-829.555 

93.004 

360.806 

47.552 

-136.604 

y shear 

859.737 

-373.360 

-297.710 

161.143 

103.938 

28.909 

zshear 

619.491 

650.008 

-132.888 

14.012 

-27.354 

72.013 

114 kn 

x shear 

-559.197 

-927.885 

106.553 

431.586 

51.905 

-138.333 

y shear 

866.258 

-389.821 

-323.750 

152.446 

110.332 

28.936 

z shear 

595.015 

983.168 

-123.824 

-49.834 

-37.713 

83.429 

128 kn 

x shear 

-665.538 

-1001.118 

150.310 

475.911 

51.352 

-122.022 

y shear 

860.691 

-422.950 

-422.327 

154.754 

105.556 

25.445 

zshear 

527.788 

1377.519 

-77.821 

-224.695 

-30.798 

36.353 

142 kn 

x shear 

-860.880 

-1019.420 

244.469 

466.649 

35.931 

-95.925 

y shear 

865.340 

-500.791 

-413.128 

225.223 

84.294 

10.720 

zshear 

522.532 

1834.928 

3.573 

-341.358 

-2.734 

4.944 
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Conversion of C81 Hub Loads For NASTRAN Input 


The sine and cosine components of the harmonics of hub shears calculated by C81 must be converted to the 
proper sign for NASTRAN. The coordinate systems for C81 and NASTRAN are shown in the figure. 

The sign change is accounted for in the phase angle representation of NASTRAN dynamic loading. Care must 
be used by the analyst when determining phase input because the phase angle equation has limitations. A ± 
180° change may be required to the predicted phase angle because of these limitations. The limitations 
stem from the inability of the arctangent function to distinguish between quadrants I and III and 
quadrants II and IV for phase angle predictions (i.e., a positive cosine, positive sine input and a - 
cosine, -sine input yield identical phase angle results when in reality they are 180° out of phase). 
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CONVERSION OF C81 HUB LOADS FOR NASTRAN INPUT 


+ Y 

RIGHT 


+ X 

FORWARD 



♦ Z 

DOWN 


+ Z UP 


+ Y RIGHT 



♦ X AFT 


C-81 COORDINATE SYSTEMS 

C-81 HARMONIC SERIES LOAD OUTPUT 


where 


/U) - a 0 + 


N 


I 


a k cos ( 2nk o> t) + sin (2nk o> t) 


NASTRAN 


fit) = TIME HISTORY DYNAMIC LOAD VECTOR 


a ,a L ,b L = INTEGRATION CONSTANTS 

O k k 

Z = HARMONIC FREQUENCY OF INTEREST (2p,4p,6p) 


NASTRAN DYNAMIC LOAD INPUT 


t = TIME 


V 2 . ,2 

X = a k + b k 

* = lan ~{z) 
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Trim to Cyclic Hub Loads (NASTRAN Format) 


ThP huh load nredictions for the 2, 4, and 6-per-rev main rotor harmonics, when converted to NASTRAN 
amplitude 1 and p^ase dynlic Toad format lor the trim-to-cycllc ("wind tunnel") cate, are shown below. 
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TRIM TO CYCLIC HUB LOADS (NASTRAN FORMAT) 


Airspeed and 
Direction 

Amplitude 

2p 

Phase 

Amplitude 

4p 

Phase 

Amplitude 

6p 

Phase 

67 kn 
x shear 
y shear 
z shear 

-908.25 

887.45 

550.40 

184.82 

103.75 

128.72 

-402.28 

370.88 

-259.77 

-18.07 

250.22 

144.50 

-94.29 

83.57 

-96.52 

143.25 

53.85 

109.19 

85 kn 
x shear 
y shear 
z shear 

-1015.63 

945.31 

-650.88 

196.91 

109.22 

85.07 

-484.72 

400.88 

-27.18 

-2.79 

-89.27 

116.70 

-156.96 

130.62 

-36.64 

161.97 

72.31 

80.61 

101 kn 
x shear 
y shear 
z shear 

-1145.42 

1024.43 

-873.07 

192.70 

105.67 

57.31 

-563.01 

448.74 

-119.46 

-11.07 

265.19 

-60.73 

-203.01 

167.44 

-39.68 

153.69 

63.90 

8.89 

114 kn 
x shear 
y shear 
z shear 

-1347.38 

1144.03 

-1104.63 

192.17 

103.48 

46.23 

-668.98 

555.92 

-140.45 

-10.34 

264.07 

-86.22 

-233.12 

190.04 

-67.09 

150.89 

62.04 

-3.64 

128 kn 
x shear 
y shear 
z shear 

-1496.91 

1200.43 

-1356.54 

184.26 

94.90 

37.57 

-763.66 

654.90 

-219.21 

-17.84 

254.54 

232.85 

-257.49 

211.01 

-19.65 

140.87 

51.69 

-31.98 

142 kn 
x shear 
yshear 
z shear 

-1897.28 

1360.04 

-1672.52 

223.20 

128.46 

27.22 

-922.69 

788.93 

-206.26 

27.02 

-63.50 

208.15 

-242.80 

210.04 

-31.73 

178.69 

95.39 

231.19 
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Full Aircraft Trim Hub Loads (NASTRAN Format) 


The hub load predictions for the 2, 4, and 6-per-rev main rotor harmonics, when converted to NASTRAN 
amplitude and phase dynamic load format for the full aircraft trim case, are shown below. 


86 



FULL AIRCRAFT TRIM HUB LOADS (NASTRAN FORMAT) 


Airspeed and 
Direction 

Amplitude 

2p 

Phase 

Amplitude 

4p 

Phase 

Amplitude 

6p 

.Phase 

67 kn 
x shear 
y shear 
z shear 

-731.22 

817.47 

-708.46 

217.95 

116.15 

123.92 

-260.38 

203.12 

-338.93 

41.56 

-50.76 

130.16 

-63.50 

50.74 

-123.50 

213.44 

128.73 

61.03 

85 kn 
x shear 
y shear 
z shear 

-848.41 

873.23 

-638.40 

213.43 

113.64 

72.40 

-309.22 

281.70 

-48.25 

23.73 

-56.82 

172.67 

-102.27 

74.47 

-57.28 

170.04 

85.02 

15.00 

101 kn 
x shear 
y shear 
z shear 

-968.56 

937.31 

-897.93 

211.08 

113.47 

43.62 

-372.60 

338.52 

-133.62 

14.45 

-61.57 

-83.98 

-144.64 

107.88 

-77.03 

160.81 

74.46 

-20.80 

114 kn 
x shear 
y shear 
z shear 

-1083.36 

949.93 

-1149.20 

211.08 

114.23 

31.18 

-444.54 

403.64 

-133.48 

13.87 

-67.81 

248.08 

-147.75 

114.06 

-91.56 

159.43 

75.30 

-24.32 

128 kn 
x shear 
y shear 
z shear 

-1202.16 

959.00 

-1475.17 

213.62 

116.17 

20.96 

-499.08 

449.79 

-237.79 

17.53 

-69.88 

199.10 

-132.39 

108.58 

-47.65 

157.18 

76.45 

-40.27 

142 kn 
x shear 
y shear 
z shear 

-1334.29 

999.80 

-1907.88 

220.18 

120.06 

15.90 

-526.81 

470.53 

-341.38 

27.65 

-61.40 

179.40 

-102.43 

84.97 

-5.65 

159.47 

82.75 

151.06 
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NASTRAN LOAD CONDITIONS 


Reaardless of the sophistication of the FEM used for the airframe vibration predictions, the ^curacy of 
the predicted response depends largely on the accuracy of the predicted rotor-induced loads transmitted t 
££& ™s“oX environment is very complex Two “r1m C1 “ n add it on 

through the control actuators and aerodynamic forces acting directly on the airframe also affect 
dynamic response. These load conditions are depicted in the figure. 
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NASTRAN LOAD CONDITIONS 


2, 4, AND 6 /REV HUB LOADS^ 

(C81)* / 

^ St, 




2 /REV CONTROL LOADS 
(OLS DATA) 


2 /REV FIN LOAD 

AT 142 KT 

(KAMAN DATA) 








*C81 HUB SHEAR PREDICTIONS 

1. TRIM TO CYCLIC ("WIND TUNNEL" CONDITION) 

2. FULL AIRCRAFT TRIM 



2P CONTROL LOADS FROM OLS DATA 


Measured 2 d boost cylinder axial loads were obtained from the OLS report (Reference 1) for application to 
the S NASTRAN fusel age ^odel f The effect of these loads on the vibratory response calculations will be 


assessed. 
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2p CONTROL LOADS FROM OLS DATA 


AIRSPEED 

F/A CYCLIC 

LATERAL CYCLIC 

COLLE 

CTIVE 

(kt) 

AMP 

PHASE 

AMP 

PHASE 

AMP 

PHASE 

67 

336.8 

139.9 

502.2 

-160.4 

294.6 

15.8 

85 

347.4 

137.2 

570.1 

-155.6 

371.6 

18.7 

101 

416.8 

122.1 

600.1 

-164.1 

557.6 

0.0 

114 

455.6 

120.5 

640.4 

-176.8 

782.9 

-1.7 

128 

560.4 

116.6 

728.3 

168.6 

925.9 

t 

P 

k> 

142 

721.5 

126.3 

815.8 

177.1 

972.8 

10.9 
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KAC Test Data for Force Determination 


2p hub shears and tail fin lateral loading at 142 kt-level flight from the Kaman AH-1G force determination 
tests (Reference 7, pp. 127 and 134) are shown below. Only the lateral fin load was applied to the 2p 
lateral response case to assess the effects on vibratory response calculations. 
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KAC TEST DATA FOR FORCE 
DETERMINATION 


FORCE 

DIRECTION 

FORCE 

DETERMINATION 

MAG 

(lb) 

PHASE 

(deg) 

VERTICAL AT HUB* 

1342. 

65 

LONGITUDINAL AT HUB* 

309. 

112 

LATERAL AT HUB* 

205. 

240 

LATERAL AT TAIL ROTOR GEARBOX 

146. 

218 


*C81 calculated hub loads were used rather than force 
determination hub loads 


93 





















j 

I 


i 

i 



in 

Ch 



PRECEDING PAGE BLANK NOT FILMED 


m kJlJL _w miuNALu aua* 


2, 4, and 6-per-Rev Flight Vibration Comparisons 

Vibration response comparisons of the coupled rotor/fuselage analyses and OLS test measurements 
are presented in this section. The comparisons are presented in the following sequence: 

1. 2, 4, and 6-per-rev vertical and lateral response comparisons (hub shears only) 

a. C81 Hub Shears - Trim-to-cyclic 

b. C81 Hub Shears - Full Trim 

2. 2-per-rev vertical and lateral response comparisons (combined loads) 

a. C81 Hub Shears (Trim-to-cyclic) + Control Loads 

b. C81 Hub Shears (Trim-to-cyclic) + Control Loads + Fin Loads 
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2, 4, AND 6-PER-REV FLIGHT VIBRATION COMPARISONS 


2, 4, AND 6 /REV VERTICAL AND LATERAL RESPONSE COMPARISONS - 
(HUB SHEARS ONLY): 

C81 HUB SHEARS - TRIM-TO-CYCLIC 
C81 HUB SHEARS - FULL TRIM 
2 /REV VERTICAL AND LATERAL RESPONSE COMPARISONS 

COMBINED 
LOADS 


C81 HUB SHEARS (TRIM-TO-CYCLIC) 


CONTROL LOADS 


LATERAL FIN LOAD AT 142 KT 
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TWO-PER-REV VERTICAL RESPONSE - HUB SHEARS ONLY 


Response calculations for six airspeeds from 67-142 kn with hub shears applied for the trim-to- 
cyclic and full aircraft trim dynamic analyses are compared for the following response points: 




FS 

1 . 

Nose 

46 

2. 

Gunner 

93 

3. 

Pilot 

148 

4. 

Engine Deck 

250 


Test data from the Operational Loads Survey Report (Reference 1) are used for comparison with the 
calculated vibrations. 

There is generally good agreement between analysis and test. 


r~D 
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AOCEL <g) ACCEL (g) 


TWO-PER-REV VERTICAL RESPONSE - HUB SHEARS ONLY 




• TEST 

TRIM TO CYCLIC 

FULL AIRCRAFT TRIM 
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TWO-PER-REV VERTICAL RESPONSE - HUB SHEARS ONLY 


Response calculations for six airspeeds from 67-142 kn with hub shears applied for the trim-to- 
cyclic and full aircraft trim dynamic analyses are compared for the following response points: 


1 . 

Tailboom Junction 

FS 

300 

2. 

Elevator 

400 

3. 

Tailboom, Aft 

485 

4. 

Tail Rotor Gearbox 

518 


Test data from the Operational Loads Survey Report (Reference 1) are used for comparison with the 
calculated vibrations. 
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TWO-PER-REV LATERAL RESPONSE - HUB SHEARS ONLY 


Response calculations for six airspeeds from 67-142 kn with hub shears applied for the trim-to- 
cyclic and full aircraft trim dynamic analyses are compared for the following response points: 




FS 

1 . 

Nose 

46 

2. 

Gunner 

93 

3. 

Pilot 

148 

4. 

Engine Deck 

250 


Test data from the Operational Loads Survey Report (Reference 1) are used for comparison with the 
calculated vibrations. 

The correlation between test and both analysis cases is poor. The calculations are much lower 
than test. Neither case predicts any significant vibration levels in the lateral direction. 
This is suspected to be due to not accounting for lateral 2/rev fin loading due to main rotor 
downwash which is evaluated in this section of the report. 
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TWO-PER-REV LATERAL RESPONSE - HUB SHEARS ONLY 



AIRSPEED (knots) AIRSPEB) (knots) 


• TEST 

TRIM TO CYCLIC 

FULL AIRCRAFT TRIM 
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TWO-PER-REV LATERAL RESPONSE - HUB SHEARS ONLY 


Response calculations for six airspeeds from 67-142 kn with hub shears applied for the trim-to- 
cyclic and full aircraft trim dynamic analyses are compared for the following response points: 




FS 

1 . 

Tail boom Junction 

300 

2. 

Elevator 

400 

3. 

Tail Fin 

521 


Test data from the Operational Loads Survey Report (Reference 1) are used for comparison with the 
calculated vibrations. 

This is again thought to be due to not accounting for the 2/rev lateral tail fin loading. 
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TWO-PER-REV LATERAL RESPONSE - HUB SHEARS ONLY 



AIRSPEED (knot*) AIRSPEB) (knot*) 


• TEST 

TRIM TO CYCLIC 

FULL AIRCRAFT TRIM 



«5 IS 105 125 145 
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TWO-PER-REV MAIN ROTOR HUB RESPONSE - HUB SHEARS ONLY 


Response calculations for six airspeeds from 67-142 kn are compared at the main rotor hub. Only 
hub shears are applied to the NASTRAN model for the trim-to-cyclic and full aircraft trim dynamic 
analyses. Test data from the Operational Loads Survey Report (Reference 1) are used for 
comparison with the calculated vibrations. 
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FOUR-PER-REV VERTICAL RESPONSE - HUB SHEARS ONLY 


Rp^nonse calculations for six airspeeds from 67-142 kn with hub shears applied for the trim-to- 
Cyclic ®nd f“l aircraft trim dynamic analyses are compared for the following response points: 




FS 

1 . 

Nose 

46 

2. 

Gunner 

93 

3. 

Pilot 

148 

4. 

Engine Deck 

250 


Test data from the Operational Loads Survey Report (Reference 1) are used for comparison with the 
calculated vibrations. 
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FOUR-PER-REV VERTICAL RESPONSE - HUB SHEARS ONLY 




• TEST 

TRIM TO CYCLIC 

FULL AIRCRAFT TRIM 









FOUR-PER-REV VERTICAL RESPONSE - HUB SHEARS ONLY 


Response calculations for six airspeeds from 67-142 kn with hub shears applied for the trim-to- 
cyclic and full aircraft trim dynamic analyses are compared for the following response points: 


FS 


1 . 

Tail boom Junction 

300 

2. 

Elevator 

400 

3. 

Tail boom. Aft 

485 

4. 

Tail Rotor Gearbox 

518 


Test data from the Operational Loads Survey Report (Reference 1) are used for comparison with the 
calculated vibrations. 
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FOUR-PER-REV LATERAL RESPONSE - HUB SHEARS ONLY 






FS 

1 . 

Nose 

46 

2. 

Gunner 

93 

3. 

Pilot 

148 

4. 

Engine Deck 

250 


Test data from the Operational Loads Survey Report (Reference 1) are used for comparison with the 
calculated vibrations. 
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FOUR-PER-REV LATERAL RESPONSE - HUB SHEARS ONLY 
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• TEST 
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FOUR-PER-REV LATERAL RESPONSE - HUB SHEARS ONLY 


Response calculations for six airspeeds from 67-142 kn with hub shears applied for the trim-to- 
cyclic and full aircraft trim dynamic analyses are compared for the following response points. 


1 . 

Tail boom Junction 

FS 

300 

2. 

Elevator 

400 

3. 

Tail Fin 

525 


Test data from the Operational Loads Survey Report (Reference 1) are used for comparison with the 
calculated vibrations. 
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FOUR-PER-REV LATERAL RESPONSE - HUB SHEARS ONLY 




• TEST 

TRIM TO CYCLIC 

FULL AIRCRAFT TRIM 
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FOUR-PER-REV MAIN ROTOR HUB RESPONSE - HUB SHEARS ONLY 


Response calculations for six airspeeds from 67-142 kn are compared at the hub. Only hub shears 
are applied to the NASTRAN model for the trim-to-cyclic and full aircraft trim dynamic analyses. 
Test data from the Operational Loads Survey Report (Reference 1) are used for comparison with the 
calculated vibrations. Note the differences in scale between the fore-aft/lateral response plots 
and the vertical response plot for proper comparison. 
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• TEST 

TRIM TO CYCLIC 

FULL AIRCRAFT TRIM 



•s 105 1*5 

AIRSPEH) (knots) 





SIX-PER-REV VERTICAL RESPONSE - HUB SHEARS ONLY 


Response calculations for six airspeeds from 67-142 kn with hub shears applied for the trim-to- 
cyclic and full aircraft trim dynamic analyses are compared for the following response points: 




FS 

1 . 

Nose 

46 

2. 

Gunner 

93 

3. 

Pilot 

148 

4. 

Engine Oeck 

250 


Test data from the Operational Loads Survey Report (Reference 1) are used for comparison with the 
calculated vibrations. 
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SIX-PER-REV VERTICAL RESPONSE - HUB SHEARS ONLY 
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• TEST 
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FULL AIRCRAFT TRIM 
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SIX-PER-REV VERTICAL RESPONSE - HUB SHEARS ONLY 


Response calculations for six airspeeds from 67-142 kn with hub shears applied forthetrim-to- 
cyclic and full aircraft trim dynamic analyses are compared for the following response points. 


1 . 

Tailboom Junction 

FS 

300 

2. 

Elevator 

400 

3. 

Tailboom, Aft 

485 

4. 

Tail Rotor Gearbox 

518 


Test data from the Operational Loads Survey Report (Reference 1) are used for comparison with the 
calculated vibrations. 
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SIX-PER-REV VERTICAL RESPONSE - HUB SHEARS ONLY 
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SIX-PER-REV LATERAL RESPONSE - HUB SHEARS ONLY 


Response calculations for six airspeeds from 67-142 kn with hub shears applied for the trim-to- 
cyclic and full aircraft trim dynamic analyses are compared for the following response points: 




FS 

1 . 

Nose 

46 

2. 

Gunner 

93 

3. 

Pilot 

148 

4. 

Engine Deck 

250 


Test data from the Operational Loads Survey Report (Reference 1) are used for comparison with the 
calculated vibrations. 
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SIX-PER-REV LATERAL RESPONSE - HUB SHEARS ONLY 


Response calculations for six airspeeds from 67-142 kn with hub shears applied for the trim-to- 
cyclic and full aircraft trim dynamic analyses are compared for the following response points: 

1. Tailboom Junction 

2. Elevator 

3. Vertical Tail Fin 

Test data from the Operational Loads Survey Report (Reference 1) are used for comparison with the 
calculated vibrations. 
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SIX-PER-REV LATERAL RESPONSE - HUB SHEARS ONLY 



TRIM TO CYCLIC 

- FULL AIRCRAFT TRIM 



185 125 145 
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SIX-PER-REV MAIN ROTOR HUB RESPONSE - HUB SHEARS ONLY 


Response calculations for six airspeeds from 67-142 kn are compared at the main rotor hub. Only 
hub shears are applied. to the NASTRAN model for the trim-to-cyclic and full aircraft trim dynamic 
analyses. Test data from the Operational Loads Survey Report (Reference 1) are used for 
comparison with the calculated vibrations. Note the differences in scale between the fore- 
aft/lateral response plots and the vertical response plot for proper comparison. 
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SIX-PER-REV MAIN ROTOR HUB RESPONSE - HUB SHEARS ONLY 
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TWO-PER-REV VERTICAL RESPONSE - HUB AND CONTROL LOADS 


Response calculations for six airspeeds from 67-142 kn with hub shears and boost cylinder loads 
applied simultaneously for the trim-to-cyclic and full aircraft trim dynamic analyses are 
compared for the following response points: 




FS 

1 . 

Nose 

46 

2. 

Gunner 

93 

3. 

Pilot 

148 

4. 

Engine Deck 

250 


Test data from the Operational Loads Survey Report (Reference 1) are used for comparison with the 
calculated vibrations. 
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TWO-PER-REV VERTICAL RESPONSE - HUB AND CONTROL LOADS 



• TEST 

TRIM TO CYCLIC 


FULL AIRCRAFT TRIM 







TWO-PER-REV VERTICAL RESPONSE - HUB AND CONTROL LOADS 


Response calculations for six airspeeds from 67-142 kn with hub shears and boost cylinder loads 
applied simultaneously for the trim-to-cyclic and full aircraft trim dynamic analyses are 
compared for the following response points: 


1 . 

Tailboom Junction 

FS 

300 

2. 

Elevator 

400 

3. 

Tailboom, Aft 

485 

4. 

Tail Rotor Gearbox 

518 


Test data from the Operational Loads Survey Report (Reference 1) are used for comparison with the 
calculated vibrations. 
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TWO-PER-REV LATERAL RESPONSE - HUB, CONTROL, AND FIN LOAD 


In addition to the hub shears and boost cylinder loads being applied simultaneously, a tail fin 
lateral load was added at 142 kn (see Page 92) for the trim-to-cyclic and full aircraft trim 
dynamic analyses. Response calculations for six airspeeds from 67-142 kn are compared for the 
following response points: 




FS 

1 . 

Nose 

46 

2. 

Gunner 

93 

3. 

Pilot 

148 

4. 

Engine Deck 

250 


Test data from the Operational Loads Survey Report (Reference 1) are used for comparison with the 
calculated vibrations. 
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TWO-PER-REV LATERAL RESPONSE - HUB, CONTROL, AND FIN LOAD 


In addition to the hub shears and boost cylinder loads being applied simultaneously, a tail fin 
lateral load was added at 142 kn (see Page 92) for the trim-to-cyclic and full aircraft trim 
dynamic analyses. Response calculations for six airspeeds from 67-142 kn are compared for the 
following response points: 




FS 

1 . 

Tail boom Junction 

300 

2. 

Elevator 

400 

3. 

Vertical Tail Fin 

521 


Test data from the Operational Loads Survey Report (Reference 1) are used for comparison with the 
calculated vibrations. 
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TWO-PER-REV MAIN ROTOR HUS RESPONSE - HUB AND CONTROL LOADS 


Response calculations for six airspeeds from 67-142 kn are compared for the hub response point 
with hub shears and boost cylinder loads applied simultaneously for the trim-to-cyclic and full 
aircraft trim dynamic analyses. Test data from the Operational Loads Survey Report (Reference 1) 
are used for comparison with the calculated vibrations. 
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Test data from an AH-1G Operational Load Survey (OLS) were used for correlation of a coupled 
rotor/fuselage analysis. Analytic predictions of hub shears (C81), OLS measured control loads, 
and a tail rotor gearbox lateral force were used to excite the NASTRAN model. The control loads 
and gearbox lateral force only apply to 2p main rotor harmonic excitation while the hub shears 
were determined for 2p, 4p, and 6p harmonics. 

The correlation was based on comparing vibration amplitudes in the lateral and vertical 
directions at selected fuselage locations for six airspeeds from 67-142 kn. 

Conclusions from the rotor/fuselage coupling analysis and flight vibration correlation study are 
as follows: 

1. An existing analysis method for coupling rotor and fuselage dynamic analyses using the BHTI 
Rotorcraft Flight Simulation Program C81 and NASTRAN demonstrated the capability to produce 
reasonable response predictions for rotorcraft. 

2. There is good agreement between calculated and measured vertical two-per -rev vibration. 
This is the predominant excitation frequency for a two-bladed teetering rotor. The 
inclusion of 2p control loads significantly affects the response level of certain response 
locations. 

3. Lateral two-per-rev vibration levels predicted by the coupled rotor/fuselage analyses were 

much lower than the measured vibrations at most of the correlation points. Main rotor 

downwash on the fin is suspected to be a factor since the correlation was noted to improve 
significantly when the two-per-rev fin load was included in the calculations for the 142 kn 
condition. 

4. Calculated and measured four- and six-per-rev vibration responses agree fairly well. It is 
not surmised that the accuracy of the analysis at these frequencies can be judged, however, 
since the vibration response at these frequencies was not strongly influenced by modes in 
close proximity to the forcing frequency. From the results of vibration test correlations, 
the airframe vibration prediction was quantitatively accurate through four-per-rev, but 
deviated from measured results significantly at six-per-rev. In addition, when exciting at 
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CONCLUSIONS (Continued) 


the main rotor hub through the pylon, the correlation obtained was poorer than that obtained 
when exciting directly on the airframe. Considering these factors, the agreement between 
test and analyses for the four- and six-per-rev responses may have been coincidental. More 
information is needed in order to judge the prediction of airframe vibration at these 
frequencies. 

Recommendations for further investigations are as follows: 

1. Vibration prediction in the four-per-rev frequency range and above needs further 
investigation. 

2. Investigate the effect of pylon dynamics on airframe vibration by a combined analytical and 
test correlation program. 

3. Investigate the main rotor two-per-rev downwash environment on the AH-1G fin. 

4. A convenient method for measuring hub shears should be developed for direct correlation with 
the analysis. 
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CONCLUSIONS 


• SUCCESSFUL DEMONSTRATION OF ROTOR/AIRFRAME COUPLING ANALYSIS 

- ROTOR ANALYSIS - C81 

- FUSELAGE ANALYSIS - NASTRAN 

• TWO-/ FOUR-, AND SIX-PER-REV CORRELATION 

- TWO-PER-REV VERTICAL GOOD 

- TWO-PER-REV LATERAL POOR (FIN DOWNWASH LOADING 
SUSPECTED) 

- FOUR- AND SIX-PER-REV FAIR 

• FURTHER WORK NEEDED 

- INVESTIGATE HIGHER FREQUENCY (>. FOUR-PER-REV) VIBRATION 
PREDICTIONS 

- INVESTIGATE PYLON DYNAMICS 

- MEASURE MAIN ROTOR DOWNWASH ON FIN 

- DEVELOP HUB SHEAR MEASUREMENT TECHNOLOGY 
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